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ABSTRACT 


The evolution of spacecraft of optimum design and ac- 
ceptable reliability requires knowledge of in-flight dynamic 
behavior during configuration changes. Systems having 
radially telescoping or folded appendages are analyzed. 
Neither system is completely amenable to closed form solu- 
tion, but analog or digital computer techniques can be ap- 
plied to specific cases. Digital techniques compare favor- 
ably with analog methods for this application. Coriolis 
effects are significant for rapid deployments. Behavior 
parameters can be related to system geometry in closed 
form but time histories require iteration. Telescoping 
systems are basically irreversible and folding systems 
cyclic but both are terminated in practice by latching de- 
vices. Solutions of the two basic systems may be rationally 
extrapolated to evaluate adequately more complex systems. 
Finally, pertinent test experience is summarized. When 
centrifugal appendage deployment is concurrent with despin, 
environmental simulation presents major problems includ- 
ing incorrect ambient gravity. Test methods are described 
with particular attention to gravity compensating techniques. 
Test results correlate acceptably with analytical predic- 
tions; however, it is concluded that analysis and testing 
should be complementary rather than mutually exclusive. 
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SIMULATION OF DEPLOYMENT DYNAMICS 
OF SPINNING SPACECRAFT 


by 

William E. Lang and George H. Honeycutt 
Goddard Space Flight Center 


INTRODUCTION 

Spin stabilized spacecraft often use the rotational energy of a constrained launch configuration 
to accomplish change to an orbital configuration of increased inertia and reduced spin rate by 
sequential or concurrent centrifugal deployment of various appendages. This paper is concerned 
with evaluation of the dynamic effects of such deployments, particularly the nonlinear transient 
phenomena generated during periods of variable inertia. This evaluation requires mathematical 
simulation of system behavior and may also require laboratory tests of system performance with 
appropriate environmental control including the negation of gravity effects. 

There has been prior recognition of appendage deployments as phenomena requiring both 
analytical and empirical study. Forsythe (References 1 and 2) and Gluck and Gale (References 
3 and 4) have reported both theoretical analysis and testing of satellites with deployable appendages. 
Either analysis or testing may, for an actual spacecraft design, involve numerous complicating 
factors which must be considered in the overall evaluation. Such factors might include complex 
mechanical linkages, distributed mass effects, initial N precession of the system, devices to assist 
or retard deployment, detent features to control or delaj deployment, frictional effects, and the 
need to consider possible failure mode performance. A more general factor, less tangible but of 
obvious practical significance, is that for any specific spacecraft design the primary objective is 
to establish confidence that it will perform acceptably or that it will or may not, with deficiencies 
adequately identified to facilitate design improvements. Also, it is common to supplement even the 
most rigorous and complete analysis w,ith some operational testing— primarily to verify quality 
control of components and assembly rather than to verify analysis. As either rigorous analyses or 
extensive test programs tend to produce much data which, however interesting, are irrelevant to 
the primary objective, along with a few values of critical significance, the purpose of this study is 
to develop methods to yield the critically significant values as quickly as possible. 

The concomitant compromise of accuracy has to be considered but usually will be tolerable 
because design practice necessarily includes an operational safety factor — a marginal design 
would normally be considered unacceptable. To be more snecific, failure modes for a satellite 
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appendage deployment system would generally be either incomplete deployment or structural 
failure during deployment or at lockup into deployed position. The precise kinematics of deploy- 
ment are often irrelevant to evaluation except insofar as they impose loads and stresses on the 
structure. 

The following presentation will first discuss general principles involving formulations based 
on energy and momentum conservation, applicable to any case of centrifugal deployment of ap- 
pendage from a spinning free body. Specific formulation defining transient behavior will then be 
developed for two basic systems. The first postulates symmetrical point masses receding radially 
from a spinning central body, and the second considers point masses constrained to circular arcs 
coplanar with the spin axis. These basic systems characterize spacecraft having radially tele- 
scoping and hinged appendages, respectively. 

Various methods can be used to solve these basic systems, and the nature of the solutions will 
be discussed. Then the applicability of these solutions to the evaluation of more complex real 
systems will be considered. It is contended that the solutions of the basic systems may suffice 
by rational extrapolation to yield adequate evaluation of many more complex systems from the 
viewpoint of critical design criteria. 

The evolution of testing facilities for subjecting spinning spacecraft to appendage deployment 
will be described and pertinent test experience will be reported. Particular attention will be given 
to the techniques which may be used to nullify gravity effects and to the significance of this factor. 
A discussion of testing facilities is included because theoretical analysis and test operations are 
considered as complementary rather than alternative evaluation techniques. 


GENERAL PRINCIPLES 

It is axiomatic that any spinning spacecraft whether of constant or varying inertia must main- 
tain constant angular momentum. The total energy of the system must also be invariant and may 
or may not include energy contributed by springs or similar devices to assist deployment or ex- 
pended in friction or deployment retarding devices. In many cases the energy involved is pre- 
ponderantly kinetic, and assumption that only kinetic energy is involved greatly simplifies analysis. 

Momentum and kinetic energy conservation yield two basic equations: 


L 


Mo = 


I T 0J 


(Angular momentum) 


( 1 ) 


2T 


Io*o 2 


I t ^t 2 




(Kinetic energy) 


( 2 ) 


The subscripts designate initial and subsequent times; I is the polar moment of inertia of the system, 
and 6 its angular velocity; "m" designates masses of parts of the system: and r and x are velocity 
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components of these masses, r being radial from the spin axis and x parallel to the spin axis. 

The third (tangential) component of absolute velocity is implicit in the I O 0 O 2 and I T 0 T 2 terms — 
therefore, Equation 2 accounts for all the systems kinetic energy. 

The specific formulation of the two basic equations will depend on the system geometry and 
will include variables &\ r and its first time derivative r. Here, r will be the expanding radius 
of gyration of the deploying part of the system or possibly some angular parameter implicitly 
defining the mechanical restraints of the system geometry during deployment. It will include a 
constant (nondeploying) inertia plus mr 2 terms for deploying masses. Velocities parallel to the 
spin axis (i.e., x components) will be expressible without introducing any other variables by 
applying a third basic principle, namely that the linear momentum of the system in any direction 
and specifically along the spin axis, is constant and may be considered zero. (Any initial linear 
velocity of the entire system is irrelevant.) Therefore 

I> = ° . (3) 

It is also true that ^ mr =0. In most cases this is intuitively obvious because of symmetry of 
the system about the spin axis (i.e., Appendages are usually designed as diametrically opposed 
similar pairs). 

While a complete and general analytical solution of the problem thus formulated might be of 
academic interest, the questions of most practical significance are: 

1. What will the spin rate be after deployment? This is obtained by substituting the deployed 
inertia value for I T in Equation 1. 

2. How much energy will be dissipated by or in the structure due to deployment? The total 

amount is equal to 1/2 (l 0 <9 0 2 or 1/2 l 0 ^o 2 ( 1_ i o /i f) where I F and d F are the final 

inertia and spin rate. The loads resulting from transfer of this energy into friction loss and/ 
or strain energy depend on the geometry and compliance of the system. The total energy is 
proportional to the square of the initial spin rate. 

3. The deployment involves transfer of momentum from the central body, which slows down, 
into the appendages, which speed up. The resulting tangential accelerations cause a rectro- 
grade "bending back” of the appendages. As this effect is not constant during the deployment, 
pertinent questions are: 

(a) For what geometric configuration is the effect maximum? 

(b) When the effect is maximum, what loads are imposed on the structure ? 

If the appendage mass "m" can be considered concentrated at a point at variable radius r 
from the spin axis, the basic Equations 1 and 2 may be written: 

= (i + mr 2 )^ 


(i + mr 0 2 ) 0 Q 


( 4 ) 
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(i + mr 0 2 ) &■ 2 


(i + mr 2 )# 2 + m (r 2 + x 2 )> 


(5) 


where I is the invariant inertia of the central body, and x will be expressible in terms of 
r , r and constant parameters of the system. 

From basic analysis of the kinematics of a particle with curvilinear motion, the tan- 
gential acceleration a t of the mass "m" is given by 

a t = r§ + 2i& (6) 


or 


a 


t 


-10 

mr 


( 7 ) 


where 20r is the so-called Coriolis acceleration. 

Equations 4 and 5 permit expressing 6 and r in terms of r and the constants. Then 
d0 /d T can be determined and also 0 since 0= (d0/dr)r. Therefore a t can be expressed as 
a function of r , and its maximum value and the value of r for which it is maximum can be 
found by solving the equation da^/dr = 0 for r and substituting the solution for a t . 

Since ma t is the tangential force acting on m and can be determined for any value of a t , 
the loads acting on the structure can be determined. The bending moment at the root 
section of the appendage will be ma^ where l is the distance from the root of the appendage 
to the point where its mass is assumed concentrated. It should be noted that if l is depend- 
ent on r , rather than constant, the maximum value of the bending moment will not coincide 
with the maximum value of a t . 

At this point the approach described has defined the angular and radial velocities and 
accelerations plus the tangential acceleration as functions of position (i.e., 0, o ', r, ?, and 
a t as functions of r). The behavior of the system is completely defined except for the time 
factor. 

4. The final question to be answered is "How rapidly does the deployment proceed and how 
long does it take?" Thus, it is necessary to express r as a time function; once this is 
done the preceding derivations of 0, f, etc., can also be converted to time functions. 

Since r = dr/dt and is a known function of r then dt/dr = l/r. Separating the variables 
and defining integration over appropriate limits, 


t 




( 8 ) 
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Therefore, the time to deploy to radius r can be determined by evaluating the integral 


J (O ' 1 dr . 


Unfortunately, even for the simple systems to be considered, this integral can be evaluated 
only by graphical integration, iterative, or computer techniques. The application of these 
general principles to two basic simple systems will now be considered. 


TYPE A SYSTEM 


This system is shown in Figure 1 and consists of a spinning spacecraft with moment of 
inertia I. It has two diametrically opposed point masses m/2 initially located at radius r 0 . At 
time t n the masses are released and thereafter slide radially out from the spin axis along 


frictionless massless guides. This system 
resembles one with telescoping appendages. 

The kinetic energy is given by 

T = -y (r 2 + r 2 0 2 ) + y l£ 2 , (9) 

which is derived from Equation 2. I is the 
moment of inertia of the spacecraft without the 
masses. The equations of motion can be ob- 
tained by using Lagrange’s equation (Reference 
5) (with no potential energy): 




Figure 1— Radially telescoping deployment system 
(type A). 


where q k is each coordinate of the system in turn. Substituting Equation 9 into Equation 10 
yields 

mr - mr <9 2 = 0 , (11) 


(i + mr 2 ) 0 + 2 mr r Q - 0 


( 12 ) 


The solution of Equations 11 and 12 for r and 8 will define the motion of the system. These equa- 
tions are nonlinear and difficult to solve in closed form. However, by reverting to the momentum 
and energy equations much useful information can be obtained relatively easily without need for a 
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complete solution. If a complete solution is desired, analog or digital computer methods must 
be used. 

Methods for partial solutions will now be presented: from energy and momentum conservation, 


(l + mr 0 2 ) 

o 

II 

HH 

+ mr 2 )0 , 

(l + mr 0 2 ) 

*„ 2 = (I 

+ mr 2 )<9 2 + mr 2 


(Note that there is no relative motion parallel to the spin axis so no k term.) 
From Equations 13 and 14 the following can be derived: 


(a) e = 



(b) 


66 

dr 


-2nr & 0 (i +mr 0 2 ) (i +mr 2 ) 


-2 


(C) r = ' (l +mr 0 2 ) 1/2 ( r2 ~ r 0 2 ) 1/2 (l + mr 2 ) ~ 1/2 , 
(d) r = t0 2 , 


(e) 0 - -2mr 6 2 (l+mr 0 2 ) 3/2 (r 2 -r 0 2 ) V2 (i +tnr 2 ) 5/2 . 


(13) 

(14) 


Using these values and Equation 6, the tangential acceleration a t can be expressed as: 

a t = r 6 + 2r 0 = 2I£ Q 2 (i + mr 2 )^ 2 (r 2 - r 0 2 ) V2 (i + mr 2 ) 5/2 


Setting da t /dr = o yields 

1+5 ror 0 2 
4m 



(15) 


and by substitution the maximum a t is 

a tmax = 0.572 I(9 0 2 (im + m 2 r 2 ) V2 . (16) 
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From a load analysis viewpoint the maximum retrograde torque developed on the central body 
(i.e.,mra t ) or the maximum bending moment (denoted as b) at the inboard end of the guide 
^i.e., ma t £r “ r 0 ^ would be more significant. Since general solutions for d(mra t )/dr = oor 
d(tna t (r - r 0 ) /dr = 0 cannot be conveniently formulated, these maxima could more easily be 
determined by graphically plotting mra t and ma t (r - r 0 ) for appropriate ranges of values for r. 

It is noteworthy that a t , mra t ', and ma t (r - r 0 ) values, which may be thought of as transient 
wrenching effects, are all basically proportional to the square of the initial spin rate. 

With phenomena of interest defined as functions of r , it remains to express r (and implicitly 
the other phenomena) as time functions. 

The time t for m to deploy to radius r is given by: 



_I_ + I 


1/2 

t = d - = ( r o e o) 1 (ji + 1 ) ( o 

rar 0 2 

\ T o) 

dr (17) 


2 

1 - 

1 


It is noteworthy that for any particular system, deployment time is inversely proportional to initial 
spin rate. It is not possible to express the integral of Equation 17 in closed form, except for the 
singular cases where either I or m are zero. However, an approximate evaluation can be made 
by graphical or iterative methods for any specific value of r 0 and the dimensionless parameter 
l/mr 0 2 . The computer solutions discussed later essentially performed this integration producing 
curves for r, f, r, d 9 0 9 etc. 

A graphical integration approach encounters the difficulty that the integrand function 



becomes infinite as r approaches r 0 . This problem can be circumvented by noting that for 
r < rj + 8 where S « r 0 , r % r 0 6> 0 2 . Therefore the time to deploy from r 0 to r 0 + 8 is essen- 
tially (2S/r 0 8 0 2 ) 1/2 and this time can be added as a constant to an integration from r = [r 0 + s] 
to the final value of r. 

If either I or m are zero in Equation 17, closed form integration is possible. For I = 0, or 
I/mr 0 2 = 0, one may derive that t = 6~ l [(r/r 0 
or I/mr 0 2 = CD , 


t = So' 1 


) 2 - lj 1/2 or r = r 0 (l + e 2 t 2 ) x/2 whereas for m = 0, 


log e 


r + Vr 2 ~ Q 
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For real systems, numerical values of I » m or m » I would be unlikely; however, all real 
systems are between these extremes as demonstrated in Figures 9 and 10 and these closed form 
expressions were used to verify generalized analog computer solutions for l/mr 0 2 = 0 and 
I / mr 0 2 = 1024 . (Computer operation precluded input of i/ mr 0 2 = <*>). For I » m or m » I the 
relationship between r and t depends only on 9 0 and r 0 . 

The behavior of any Type A system is therefore predictable with relative ease even without 
the computer techniques which will be described later. At this point, it is interesting and instructive 
to examine the nature of solutions for a specific case which are presented in Figures 2 and 3. For 
this example, I = 10.5 slug-ft 2 , m = 0.25 slugs, r 0 = 1 ft., d Q = 4.82 rad/sec. Figure 2 presents 



• I V Tr / 

Figure 2— Type A system ?, a f , 0, B and time versus r. 

Figure 3— Type A system r, a t , 6 and B versus time. 

r, a f , b y B (the maximum bending moment at the 
inboard end of the guide, which is ma t [r - r 0 ]) 

and time, t, as functions of r . r starts at r 0 increases to a maximum and then decreases. a t 
starts at zero, increases to a maximum defined by Equation 16, and then decreases, 0 decreases 
as the system inertia increases. B rises from zero to a maximum and then decreases. Figure 3 
presents the same parameters as time functions. The various maximum are nonsimultaneous and 
occur at different values of r. 

TYPE B SYSTEM 

This system is shown in Figure 4 and consists of a spinning spacecraft as before, but this time 
the masses are attached to massless booms of length l which in turn are attached at radius r 0 . 

At tim e t 0 the booms are released and thereafter swing out from the spin axis. The central body 
has mass M . The angle between the booms and the spin axis is denoted a and is initially zero. Other 
symbols retain their significance in the prior discussion. 
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The kinetic energy of this system is ex- 
pressed by 

2T = m[^ 2 a 2 + (r 0 + sin a) 2 Q 2 + 2a Z sin a x + x 2 J 

+ Mx 2 + I<9 2 f (18) 

where x = - m l a sin a/M + m and is the linear 
velocity of m. 

If M is considered to be » m , (see Ap- 
pendix A) then the kinetic energy of this system 
is simplified: 



Figure 4— Hinged deployment system (Type B). 


2T 



+ (r Q + 'i sin a] 


< 9 2 ] + ie 2 


(19) 


Using Lagrange's Equation 10 with Equation 19: 


(i + mr 2 )&" + 2raa0r-tcos a - 0 


-t 2 a - r -£ cos ad 2 - 0 . 


( 20 ) 

( 21 ) 


where 


r = r„ + •i sin a 


These equations are obviously more difficult than those for the Type A system but can still be 
solved with the help of an analog or digital computer. The analog computer used to solve the Type A 
system did not have enough multipliers for the Type B system equations so only the digital solution 
was obtained. As with the Type A system, it is possible to use basic principles to obtain partial 
solutions of considerable value. 

This system differs from Type A in two basic ways. Both M and m develop velocity components 
parallel to the spin axis; therefore, a rigorous formulation has to include the interchange of linear 
energy and momentum between them, consistent with the net energy and momentum of the whole 
system remaining constant. Also, while system A is an open ended, irreversible phenomenon (the 
sliding weights would recede without theoretical limit with the spin rate decreasing asymptotically 


9 


to zero), system B is basically periodic. The arms would flap repetitively through a 180° arc, but 
it will suffice to limit the study to the first 90° of travel. 

Here a choice of geometric variable exists: either the radius r to the mass m, or the angle a. 
Either one implicitly defines the other since r = r 0 + l sin a. 

In terms of r , the basic equations can be derived so that 


(! + mr o 2 ) 

0 Q - (i + mr 2 )£ , 

(22) 

(i + mr 0 2 ) 9 2 = (i + mr 2 ) 

m -6 2 r 2 T m / r_r o\ 2 

"<■. L > t JJ 

(23) 


These equations are based on Equations 4 and 5 with due consideration to the fact that ^mx = 0 
(Equation 3) and the geometric constraints of the system. 

In terms of a, the equivalent equations are 


(i + mr 0 2 ) 6 0 



+ 4 sin a 



(24) 


(i + mr 0 2 ) Qq - I + m(r 0 + t sin a ) 2 6 2 + m Z 2 a 2 ^1 - jj? m sin 2 aj 


(25) 


If M » m, as in many actual systems, the constant m/(M + m) might reasonably be neglected, since 
it is essentially zero; but even then the formulation is harder to work with than that for the Type A 
system. The implications of the relative magnitudes of M and m are discussed in Appendix A. 

The tangential acceleration a t can be expressed as 


a 


t 


2I^ 0 2 (1 + mr 0 2 ) 3/2 

<1 (i + mr 2 ) 5/2 



(26) 


The expression for a t in terms of a is more cumbersome. The maximum value of a t may be found 
by trial substitution of values for r, since da t /dr = 0 is not amenable to solution. For constants 
representative of many real systems, a troax occurs between a = 20° and a = 30? The retrograde 
torque on the central body is mra t and its maximum will not coincide with a tmax . However, the 
maximum bending moment across the hinge will coincide with a tmax and will equal m4a tmax . 

The wrenching effects are, as with Type A systems, basically proportional to the square of 
initial spin rate. 
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The time t for m to deploy to radius r - r t (limited by conditions that r t 
* 90°) is given by: 



It is noteworthy that, as for the Type A system, deployment time is inversely proportional to initial 
spin rate. Integration between the limits r 0 and r t = r 0 + l would give the time to deploy through 
90°. 


In terms of <x, the time to deploy to angle a t ^ 90° is 



Integration between limits of zero and a t = n/2 would give the time to deploy through 90° 


These time integrals cannot be evaluated in closed form. It is necessary to use graphical or 
iterative methods to obtain estimates for specific system parameters. In contrast to the Type A 
system, which can be characterized by a single dimensionless parameter (l/mr 0 2 ) the Type B 
system requires three such parameters, namely l/mr 0 2 , ^/r 0 , and m/M + m and they are mutually 
independent. Therefore, though a computer could 
be programmed to solve any particular system, 
a generalized solution, like that obtained for 
the Type A system by plotting families of 
curves for a wide range of l/mr 0 2 values, is 
not feasible. 

Figure 5 compares the time history of r 
for a Type A system having the constants used 
in the previous example with the time histories 
of r and a for an equivalent Type B system 
having the same constants plus values of 4 feet 
for l and a> for M. The B system solutions 
were obtained by graphical and iterative meth- 
°ds and later confirmed by a digital computer Figure 5 _ r versu$ Hme COfnpared for Type A and Type B 
per Figures 12, 13, and 14 (described in a later systems. 


90 [ 
80 h 

70 
60 
50 
40 
30 
20 
10 
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section). The B system deployment is slower but minimally so for the period immediately after 
release when a is small. 

If M is assumed zero, the equations revert to the equivalent equations for the A system; and 
a B system withM = 0 would in fact behave as an A system with the Type A history of r, A B 
system solution was also obtained and plotted for M = 10 m = 2. 5 slugs. The resulting r versus 
time curve (not shown) deviated only slightly from the B system curve shown and moved closer to 
the A system curve with a final deployment time of 0.6135 seconds rather than 0.6232 seconds for 
M = a>. The magnitude of M also affects the tangential acceleration T (per Equation 26). T is less 
for infinite M than for finite M, but the difference is very minor for M i 10m. The relative mass 
factor is discussed further in Appendix A. 


APPLICATION OF ANALOG COMPUTERS 
Type A System 

The solution of Equations 11 and 12 can readily be obtained (Reference 6) using any general 
purpose analog computer providing the computer has enough amplifiers, integrators, multipliers, 
etc, (The analog computer used for all of the problems discussed here was the Electronic Asso- 
ciates Inc. PACE TR-48.) 

Since Equations 11 and 12 reasonably describe the motion of a number of appendage deployment 
systems, it was decided to nondimensionalize the equations and get a general solution in the form 
of families of curves which could be entered directly for approximate solutions. This was ac- 
complished by introducing the following non-dimensional variables*: 


r - 0 Q t (non-dimensional time) 


R - — (non-dimensional radius) 


& 

$ 


dR 

dr 


d 2 R 

dr 2 


r o 


*o * 2 


§ = 
0 = 


de 

dr 


d 2 Q 

dr 2 


_e_ 

e' 


*o 2 


(5 = 0 ) 


t£is inherently non-dimensional, but it is desirable to use different symbols for all variables in the generalized problem with the 
* superscripts denoting derivatives with respect to non-dimensional time r. 
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Substituting these relationships into Equations 11 and 12 yields 


S = RS’ 


( 29 ) 



+ R 2 j 


8 + 2rr3 


0 


(30) 


Equations 29 and 30 consist of dimensionless variables and one constant term I/mr 0 2 . This constant 
term defines the physical parameters of the spacecraft system, i.e., I, m, and r 0 . Thus by solving 
equations 29 and 30 for a wide range of l/mr 0 2 values one can obtain a family of curves for each 



Figure 6— Analog computer schematic for generalized 
Type A system. 


non-dimensional variable which can be transposed 
into the actual variables for any spacecraft that 
resembles the Type A system. 

Figure 6 shows the programming schematic 
for the solution of Equations 29 and 30. As a 
check on the generalized solution a specific 
problem was also programmed using the same 
physical parameters used in the example of the 
preceding alternative method of solution; 

I = 10.5 slug-ft 2 , r 0 = 1 foot, m = 0.25 slugs, 
and 0 o = 4.82 rad/sec., Figures 7 and 8 show 
the analog computer solution for this specific 
problem. 

Figures 9 and 10 show the generalized 
analog computer solution for the families of R 
and @ for values of I/mr 0 2 from 0 to 1024. 




Figure 7— Specific analog computer solution for r, r, 
and r (Type A system). 


Figure 8— Specific analog computer solution for 6 , 6 , 
and 6 (Type A system). 
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Figure 9— Generalized analog computer solution for R Figure 10— Generalized analog computer solution for 0 
(Type A system). (Type A system). 


Using the generalized curves and working back to obtain the actual variables for the specific 
problem will give the same results, thus verifying the generalized solution. The results also 
correspond closely with those obtained for the same problem by the methods previously discussed, 
and shown in Figure 3. 

Type B System 

It would be interesting and informative to perform the same kind of analog computer solution 
for the Type B system. However, Equations 20 and 21 are more complicated than Equations 11 and 
12 and any attempt to non-dimensionalize them results in the several independent constants, men- 
tioned previously, which would made the number of family curves prohibitively large. In fact, even 
to solve a specific problem requires a number of sine -cosine generators and more multipliers than 
needed for the Type A system equations. The TR-48 computer did not have enough multipliers to 
solve this problem and since the MIMIC. Digital Simulation Program (Reference 7) was available 
no attempt was made to solve the equations with an analog computer. However, with a larger analog 
computer a specific spacecraft deployment problem could be readily handled. 


APPLICATION OF DIGITAL COMPUTERS 

The equations of motion for both Type A and Type B systems can be easily solved using a 
digital computer by employing the MIMIC Digital Simulation Program explained in Reference 7. 
T his program was developed at Wright- Patterson Air Force Base to facilitate the application of 
digital computers to engineering problems. MIMIC is a program whose input language endows a 
digital computer, from the viewpoint of the user, with some significant advantages of an analog 
computer, in that the input relates directly to the physical nature of the problem, while virtually 
eliminating the time and amplitude scaling needed for analog computation. 

The Type A system with the same exemplary parametric values previously used was also 
solved using the MIMIC Program and an IBM 7094 computer. The resulting answers were virtually 
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identical and the plots of the MIMIC results are therefore the same as those shown in Figures 7 and 
8. Figure 11 shows the input required by the MIMIC Program to solve the Type A system per Equa- 
tions 11 and 12. This is the necessary input form with one punched card per line. 

Figure 12 shows the input to the MIMIC Program for solution of a Type B system. Equations 
17 and 18. The results of the Type B system solution are shown in Figures 13 and 14. The input 
values were the same as in the previous examples with = 4 feet, and M assumed infinite. A finite 
value of M could be considered without excessive complication. 


The MIMIC Program also lends itself to the solution of concurrent deployment of dissimilar 
appendages. The equations of motion would be more complicated and the computer input longer 
but the solution would still be much more rapid than by other methods. If some of the appendages 
are released sooner than others or locked in place while others are still deploying, then the problem 


could be solved in steps. This means the initial conditions and equations of motion would change 
at each abrupt change caused by appendages releasing or appendages locking in place. The MIMIC 
Program is a powerful and useful tool in the 


solution of the spinning spacecraft problem. 


SPINNING BODY WITH SLIDING WEIGHTS 


CON (M, I, DT) 

2DR = R* 1DTH *1DTH 
1DR = INT(2DR,0. ) 

R = INT( 1DR, 1.) 

A = -2.*M*R/(I + M*R*R) 

2DTH = A* 1DTH * 1DR 
1DTH = INT ( 2DTH, 4.85) 

TH = INT(1DTH,0.) 

FIN(T,2. ) 

HDR( TIME, R, RDOT, RDDOT, THDOT, THDDOT ) 
HDR 

OUT(T, R, 1DR, 2DR, 1DTH, 2DTH ) 

END 


Figure 11— Digital computer program for Type A system. 



FOLDED APPENDAGES DEPLOYMENT 


CON (I, L, R, DT) 

PAR(M) 

B = R+L * SIN (A) 

2DA = B*COS(A)« IDTH * 1DTH/L 

IDA = INT(2DA,0. ) 

A = INT ( IDA, 0. ) 

2DTH = -2.* M *1DA * 1DTH* B*L*COS(A)/( l + M»B* B) 
1DTH = INT(2DTH,5. ) 

TH = INT ( 1DTH,0. ) 

FIN(T,2. ) 

HDR(TIME, ALPHA, ADOT,A2DOT, THDOT, TH2DOT) 
HDR 

OUT ( T, A, 1 DA, 2DA, 1 DTH, 2DTH ) 

END 


Figure 12— Digital computer program for Type B system. 



Figure 13— Digital computer solution for a, a, and a. 
(Type B system). 


Figure 14— Digital computer solution for 6 and 9 (Type 
B system). 
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COMPARATIVE EVALUATION 


At this point it is pertinent to note that the same specific Type A problem was solved by three 
different and independent methods - the quasi-classical, with an analog computer and with a digital 
computer. All methods yielded essentially identical results. The actual system which led to the 
choice of these specific parameters also behaved approximately as predicted during test deployment 
(see test experience). 


EVALUATION OF ACTUAL SPACECRAFT SYSTEMS 

At this point it is pertinent to consider pragmatically the value of the preceding material. 

Real systems are less simple than those considered though deploying elements are often either 
essentially free to recede radially from the spin axis (Type A system) or essentially constrained 
to rotate out about fixed axes (Type B system). From discussion of these simple systems, the 
prospect of attempting a classical solution for a real system is discouraging especially as practical 
objectives often require only approximate answers to specific questions. Although even extremely 
complicated systems could be solved iteratively by sufficiently sophisticated analysis and computer 
facilities, the effort involved might well be inappropriate if acceptably accurate solutions leading 
to the same practical conclusions could be obtained relatively easily and quickly. 

It is suggested that a sufficiently accurate and probably conservative (i.e., more severe than 
actual) estimate of dynamic behavior can be made by assuming that a real system has a functionally 
equivalent Type A or B system. Such an educated guess is certainly preferable to neglecting a 
possible problem area. It seems reasonable to postulate a rational equivalent model for many 
actual configurations. 


Examples 

1. Appendages which telescope radially outward. These are basically Type A systems and 
could be analyzed as such by assuming the mass of sliding elements concentrated at the 
mass center. 

2. Single hinged appendages folded back against final stage booster, each as solar paddle 
arrays on many spacecraft. These are basically Type B systems with distributed de- 
ploying element mass. The total mass of each appendage could be considered concentrated 
at its mass center. Then the deployment could be analyzed as a Type B system. 

3. Double hinged appendages such as shown in Figure 15. The validity of a simplified model 
for this kind of configuration is more dubious and controversial and largely depends on the 
relative mass and inertia of the inboard and outboard elements. In general, the inboard 
element is a Type B system while the outer element tends to behave like a Type A system, 
especially if its mass greatly exceeds that of the inner element. A Type A deployment is 
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inherently more violent than a Type B 
deployment and if a Type A system, with 
deploying mass equal to the total append- 
age mass, is assumed, conservative es- 
timates of angular deceleration and de- 
ployment time will be obtained (i. e., 

Angular deceleration will be greater 
and deployment time less than for the 
actual system). The results would give 
some insight as to whether the transient 
dynamic effects during deployment might 
be significantly detrimental to the sys- 
tem. Before actual tests of a system as 
shown in Figure 15, behavior was pre- 
dicted by assuming a Type A system and 
the results obtained for duration of deployment and maximum angular deceleration 
were reasonably compatible with observed behavior during the tests. 

4. Concurrent deployment of dissimilar appendages such as a set of inertia booms and a set 
of solar paddles. A system of this type requires consideration of the extent and nature of 
dissimilarity. For minor differences, a possible approach might be to assume that (1) the 
booms deploy from the initial spin rate with the paddles initially fully deployed and then 
(2) the paddles deploy from the initial spin rate with the booms initially fully deployed. 

For each case, develop time histories of accelerations, forces and torques, including the 
effects on the deployed paddles of boom deployment and vice versa. Then, assume (1) and 
(2) to be concurrent and add the two time histories for each parameter. The results would 
be conservative as the hypothetical situation involves greater energy and momentum transfer 
than the actual concurrent deployment. 

5. Systems having springs to assist deployment and/or damping devices to retard deployment. 

A preliminary step would be to compare the energy dissipated in deployment (i. e., 

E = (1/2) (i 0 0 o 2 )(l - l 0 /l F ) with the energy supplied by the spring (+e) or absorbed by the 
dampers (-e) during deployment. If "E" considerably exceeded "e", it would suffice to add 
"e" to "E" (for springs) or subtract "e" from "E" (for dampers) and calculate a new equiv- 
alent initial spin rate as [# 0 2 ± 2e/l 0 (l - I 0 /I F )] 1/2 • Then the system could be analyzed 
for the modified initial spin rate as if the springs and/or dampers were not present. 

Forsythe 2 has applied similar logic to the problem of simulating zero gravity appendage 
deployment. 

It should be emphasized that these examples are given only to illustrate the kind of logic that 
might be usefully applied with appropriate discretion. Any particular system should obviously be 
studied with full and complete knowledge and consideration of all germane factors. Also, this kind 
of "ballpark" analysis is more appropriate in conjunction with experimental testing of system 





FOLDED DURING DEPLOYMENT DEPLOYED 

Figure 15— Gravity negation for double hinged boom 
deployment. 
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behavior. If the system cannot be tested, a more rigorous analysis would be desirable though 
even then the methods discussed here might establish sufficient confidence in operational 
reliability to omit further study - or show that the design is fundamentally unreliable. 

Testing Experience 

Various spacecraft have been subjected to appendage deployment under conditions intended to 
adequately simulate orbital situations. The most intractable problem is probably the elimination 
of unrealistic external restraints, particularly those imposed by the ambient gravitational field. 
Four methods have been used to overcome this problem. 

The first and certainly the best, when it can be used, is to have the deployment occur under 
"free fall" conditions. Unfortunately it is only feasible for short test durations for obvious reasons. 
The larger and more massive systems characteristically require longer times for deployment to 
occur and should not be subjected to the high accelerations which would necessarily be imposed by 
injection into, and recovery from, extended free fall trajectory. 

The second method, of limited value but having the merit of simplicity, is based on determining 
the kinetic energy loss due to deployment which is (1/2) (l 0 $ 0 2 ) l 1 " i o/ i f) and assuming that this is 
absorbed by the structure at lockup into deployed configuration and appropriately divided between 
the various structural elements. Then each element is subjected to a shock pulse of appropriate 
energy content. This can often be done by setting up each appendage so that gravity causes its 
deployment with supplementary energy input from springs if necessary. The entire system is not 
subjected to spin, and the kinematics of deployment, including the possibly critically significant 
Coriolis effects are not simulated. However, such subsystem tests can reveal the ability of the 
structure to withstand "lock in" shock loads. 

The third method, described by Forsythe (Reference 2) is to "overspin" the system to (6 0 + h) 
so that the excess available energy which is (1/2) I 0 (2 §£„ + S 2 ) (l - l 0 /l F ), equals, and therefore 
cancels, the gravity potential energy required to deploy the appendages. It is easy to apply and 
proper in some cases but for many cases it is entirely incorrect. For example, for a system 
shown in Figure 15 (with spin axis vertical) overspinning would cause excessively violent deploy- 
ment of the outboard elements and possible failure to lock into fully deployed position. Also, for a 
basic Type A system, telescoping radially in a horizontal plane, no overspin would be indicated, 
but the weight of the deploying elements could well cause significant frictional resistance to 
deployment. Therefore this method should be used only with due consideration as to how well it 
actually duplicates orbital conditions. At best, it imposes equivalent "lock in" shock loads; the 
dynamics of deployment will always differ from those in orbit to some extent. 

The fourth method is to actually provide a "counterforce” to each deploying element to cancel 
its weight. Elastic elements are commonly used, and they may either equate elastic energy with 
appropriate potential energy with a substantial force variation during deployment or, in more 
sophisticated applications, have a low enough spring rate (normally by virtue of length) that the 
force during deployment is relatively constant and essentially vertical, thereby approximating 
cancellation of gravity during deployment. 
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Figure 15 illustrates a test arrangement used to achieve compensation of gravity to less than 
a net variation from -0.03 g to +0.03 g (by actual measurement) during deployment of double 
hinged appendages. The elastic elements used were cis-polyisoprene (natural gum rubber of high 
purity, as used in surgical tubing). This material was selected as having superior elastic proper- 
ties and was used at elongations for minimum percentage deviation from tensions equal to the 
weights of the suspended elements. The radii to the upper attachment of the suspension elements 
were determined for minimum departure from verticality of their attachment to the appendage 
elements, after allowing for their spin-induced centrifugal deflection. The lower attachments 
were located at the mass center of each appendage element. Figure 15 shows a rigid extension 
of the spacecraft (to a height about 35 ft. above it) to drive the suspension beam - in a more recent 
application there was no physical connection except the elastic elements and the suspension beam 
remained oriented with the spacecraft within acceptable deviations. 

It is believed that this test achieved a close approximation of zero gravity deployment. The 
total counterforce, as monitored by a load cell above the suspension beam, deviated only 6 percent 
during deployment. Unrealistic horizontal components of the suspension force did not, from high 
speed movie coverage, appear to exceed 7 percent of the tension in the element. 

The test results tended to confirm predictions based on the foregoing analysis. A Type A 
system was assumed as the majority of the appendage mass was in the outer elements. The movies 
showed that the outer element mass center did, in fact, tend to recede radially in a horizontal plane 
from the spin axis. The duration of the deployment transient was approximately as predicted. The 
maximum angular deceleration from monitoring accelerometers installed tangentially on the central 
body was also approximately as predicted. The appendages (fiberglass booms) showed extreme 
deflections and hinge failure occurred during an overspeed (failure mode) test. The failure occurred 
during deployment and before iatch in and was observably due to the Coriolis wrenching effect. Ac- 
celerometers were installed on the appendage elements to monitor transient accelerations; but the 
results were difficult to interpret because of the failure and because of the transient shocks induced 
by pyrotechnic boom release devices plus other high frequency indications of indeterminate origin. 

To summarize the test results (including some not cited) tend to confirm the validity and utility 
of the preceding analysis. The accuracy of either analysis or testing cannot be precisely evaluated 
because the analysis applies to simplified and idealized systems and because the test systems, 
besides differing from the idealized models, were subject to imperfect environmental simulation. 
Since the analysis is based on incontrovertible basic principles, the computational errors would 
be minor, and the environmental simulation error was demonstrably minor, discrepancies between 
predictions and test results could be considered largely due to idealizing of the actual systems. 

Since these were no major discrepancies, the idealizing assumptions were apparently essentially 
valid. 

CONCLUSIONS 

The most significant conclusion is that the deployment dynamics of spinning spacecraft can be 
adequately simulated either analytically or by testing, but for maximum confidence the approaches 
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should be complementary. Also, the solutions of simple systems, as described, can provide 
valuable insight into the dynamic behavior of more complex systems. 

Both analysis or testing are subject to some degree of inaccuracy, the analysis being subject 
to simplifying assumptions and imprecise parametric data, and testing to imperfect simulation. 
Therefore either analysis or testing should demonstrate adequacy for conditions moderately more 
severe than those anticipated in flight. If this is done then the inaccuracies are tolerable, provided 
they are recognized. 

Both analog and digital computers are valuable tools for analysis and the relative ease with 
which useful results may be deduced without recourse to computer facilities suggests that this 
approach should not be neglected even if computer facilities are available. It affords a useful 
concept of the nature of the dynamic phenomena involved. 
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Appendix A 


Influence of Central Body to 
Appendage Mass Ratio on Deployment Dynamics 

For tests involving components of motion of the deploying appendages parallel to the spin axis, 
there has been consideration of whether the support and restraint of the central body of the space- 
craft, which effectively couples it to the comparatively infinite mass of the earth, compromises 
environmental simulation appreciably since in orbit (or in free fall tests) the central body mass is 
limited and the entire system is a free body. This general topic is discussed in this appendix and 
it is concluded that this factor can reasonably be ignored for the usual case where the central body 
mass in orbit is considerably greater than the total mass of deploying appendages. 

Referring to the discussion of Type B systems, the orbital mass ratio is M/m. For Type A 
systems this ratio is immaterial, and, as previously mentioned, if M = M/m = 0, a Type B system 
behaves like a Type A system with faster deployment and higher values of tangential acceleration 
a t . From Equation 26 the highest possible a t with M/m = o exceeds the lowest possible a t with 
M/m = oo by a factor of [4 2 /4 2 - (r - r 0 ) 2 ] 1/2 . For an arbitrarily chosen partially deployed position 
where r - r 0 = 4/2 this factor would have a numerical value of 1.156. However, a similar compar- 
ison between M/m = 1 and M/m = <» gives a numerical value of 1.069 and between M/m = 10 and 
M/m = oo gives a numerical value of 1.011. The effect of mass ratio on deployment time is simi- 
larly minor for M/m > 1 (from Equation 27). Deployment is slower for a greater value of M/m. 

Therefore the support and restraint of the central body actually results in an undertest so far 
as the transient effects during deployment are concerned but to a minor degree for usual mass 
ratios (e.g., only 1.1 percent for M/m = 10). 

The other factor to be evaluated is the shock effect at lockup. From basic principles it is 
known that this effect has a constant energy content of (1/2) i 0 $ 0 2 (l - i 0 /i f ), regardless of mass 
ratio, and since it is a collision between elements which terminate deployment (though energy is 
dissipated in compliance of the entire system) the question is essentially whether the "damage 
effect" of collision between two bodies depends only on the impact energy content or whether it 
is also influenced by the relative masses of the bodies. Their absolute and relative compliance, 
resilience and surface hardness are also pertinent factors but an adequate review of all these 
factors is a complex general problem beyond the scope of this paper. However it should be 
mentioned that precise orientation of appendages after deployment is often a design objective. 

If a collision between two bodies has to dissipate a fixed amount of energy their relative 
velocity at impact is maximum when their masses are equal and minimum when one mass is 
infinite. The maximum relative velocity exceeds the minimum by a factor of /2. Between these 
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two extremes, the relative velocity for a mass ratio of 10 is much closer to that for an infinite 
mass ratio than to that for a mass ratio of unity. However, the momentum change (i.e., impulse) 
to each body is greater for high mass ratio and lower relative velocity. 

It follows that, for the appendage lock-in situation, an infinite (coupled to earth) mass ratio will 
product a shock pulse of longer duration, but lower onset rate, than a true orbital mass ratio. The 
former sees more momentum change but at a slower rate. Both develop the same maximum forces 
and therefore the same stresses if elastic limits are not exceeded. These facts, though of academic 
interest, hardly warrant any definitive conclusion, except that the "test versus orbit" comparison 
of this factor for most systems would show it as probably insignificant as neither onset rate nor 
impulse would differ appreciably. 

Therefore, restraint of the central body during test can usually be ignored for evaluation 
purposes. Also, the central body mass could usually be considered infinite for analytical 
evaluation. 
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